Challenges for Kinetic Unified Dark Matter 



in 
o 
o 

(N 

Or 
<D 
00 

(N 
(N 

(N 
> 

o : 

in . 
o ■ 

Or 

6 ■ 
a: 

C/3 . 



13 



Dimitrios Giannakis 1 and Wayne Hu 2 
1 Department of Physics, University of Chicago, Chicago IL 60637 
2 Department of Astronomy and Astrophysics, University of Chicago, Chicago IL 60637 

Given that the dark matter and dark energy in the universe affect cosmological observables only 
gravitationally, their phenomenology may be described by a single stress energy tensor. True uni- 
fication however requires a theory that reproduces the successful phenomenology of ACDM and 
that requirement places specific constraints on the stress structure of the matter. We show that a 
recently proposed unification through an offset quadratic kinetic term for a scalar field is exactly 
equivalent to a fluid with a closed-form barotropic equation of state plus cosmological constant. 
The finite pressure at high densities introduces a cutoff in the linear power spectrum, which may 
alleviate the dark matter substructure problem; we provide a convenient fitting function for such 
studies. Given that sufficient power must remain to reionize the universe, the equation of state 
today is nonrelativistic with pap 2 and a Jeans scale in the parsec regime for all relevant densities. 
Structure may then be evolved into the nonlinear regime with standard hydrodynamic techniques. 
In fact, the model is equivalent to the well-studied collisional dark matter with negligible mean free 
path. If recent observations of the triaxiality of dark matter halos and ram pressure stripping in 
galaxy clusters are confirmed, this model will be ruled out. 



I. INTRODUCTION 

The standard model for cosmology contains dark en- 
ergy in the form of a cosmological constant and cold dark 
matter (ACDM). While this model explains most obser- 
vations to date, the dark matter and the dark energy are 
only observed gravitationally. It thus remains possible 
that they are a single phenomenon, sometimes called in 
the literature unified dark matter (UDM). 

A candidate for UDM that has been extensively dis- 
cussed in the literature is the generalized Chaplygin gas. 
Here the UDM is taken to be a perfect fluid with an equa- 
tion of state p = —A/p a , with A > and < a < 1; 
models with a = 1 0,12] are particularly attractive since 
they represent a Born-Infeld scalar field [3, B that ad- 
mits a string theory interpretation (e.g. 0, @)- This 
class of equations of state leads to a background evolu- 
tion that transitions from a matter dominated phase at 
early times to a de Sitter phase in the asymptotic future 
in agreement with distance measures 0, H H ■ 

Unfortunately, in linear perturbation theory Chaply- 
gin UDM produces oscillations or exponential growth in 
the matter power spectrum unless a — > 01 1 even after 
the effect of baryons has been included. The exact 
phenomenology of these models once nonlinear structures 
have formed is complicated due to the fact that the pres- 
sure decreases with increasing energy density (see e.g. 
[l2L Il3|). It is clear, however, that these models differ 
substantially from the successful ACDM model and are 
unlikely to be viable in their simplest form. 

These problems with Chaplygin gas UDM stem from 
very general considerations and serve as a guide to build- 
ing more successful UDM models. Fundamentally, a 
UDM model that mimics the joint stress energy tensor 
of ACDM will reproduce all of its observational successes 
[Tij . Since the UDM stress-energy tensor is covariantly 
conserved, mimicry in the spatial stresses is all that need 



be explicitly enforced p"o| . 

Deviations from perfect mimicry must be small in a 
very specific sense. The basic requirement for mimicry is 
that spatial fluctuations in the UDM stresses be negligi- 
ble in comparison to those in its energy density 5p, even 
though the homogeneous or background stress must be 
of order — p to accelerate the expansion 0, 0]. More 
specifically, the ratio between the fluctuating stress and 
the energy density must be of order the square of that 
between the momentum and energy densities, or a non- 
relativistic velocity squared. This disparity implies that 
the total pressure is effectively decoupled from the total 
den sity , i.e. that the pressure is effectively non-adiabatic 
pTot flo . fl7| . For a generalized Chaplygin gas this require- 
ment is only satisfied as a — > 0. 

In the linear regime, where the momentum density it- 
self is a perturbation, fluctuations in the spatial stresses 
must then be a second order effect quantified by a van- 
ishingly small sound speed and anisotropic stress. The 
internal dynamics of dark matter halos place even more 
stringent requirements on mimicry. Here the UDM stress 
energy tensor must embed the information from multiple 
momentum streams in the dark matter, i.e. the full phase 
space structure of cold dark matter. A successful UDM 
model will be at least as complex as ACDM in order to 
reproduce its rich phenomenology detracting somewhat 
from its appeal. 

The prototype for a successful UDM model is a scalar 
field coherently oscillating in an offset quadratic potential 
|14|: it satisfies mimicry requirements all the way down 
to the de-Broglie wavelength of the field 0|. It does not, 
however, explain the relationship between the 10 -33 eV 
flatness of the offset scale and the > 1 -22 eV mass scale. 

Recently, Scherrer 0] introduced a fc-essence 20] vari- 
ant of scalar field UDM that also succeeds in linear 
theory. Here the quadratic potential is replaced by a 
quadratic kinetic term (see 21] for a possible physical 
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motivation) . An important feature of this model is that 
its stress-energy tensor is purely kinetic (see [52 for a 
generalization). We shall show that its equation of state 
is consequently both barotropic and expressible in closed 
form. We are consequently able to reduce this kinetic 
form of UDM to collisional dark matter with a negligible 
mean free path. Non-linear structure formation for this 
case has been well-studied in the literature |23t l24j . 

The outline of the paper is as follows. We begin in 
'HTl with a brief review of the basic properties of scalar 
field UDM. We derive an exact, fully non-linear, equation 
of state for the purely kinetic UDM model [n| in AHII 
We then calculate the evolution of linear density fluctu- 
ations in 9IVI We place constraints on the model from 
the suppression of the low mass dark matter halo abun- 
dance in Sj^and discuss opportunities and challenges for 
the model from halo substructure in WII We conclude 
in WTfl 

Throughout we will compare the UDM model to a 
ACDM model with a total non-relativistic matter den- 



sity Q m h 2 



Q r h 2 



n h h 2 



0.14, baryon density 



il b h 2 — 0.024, in a flat cosmology with Q\ = 0.73 and 
scale invariant adiabatic initial conditions. 



II. SCALAR FIELD UDM 

In this section, we briefly review the relationship be- 
tween scalar fields and perfect fluids and show why scalar 
fields whose potential and kinetic terms possess a non- 
vanishing minimum value satisfy the basic requirements 
of a successful UDM model. The minimum acts as a de- 
coupled source of negative pressure with an equation of 
state p/p = — 1 for the background that does not con- 
tribute to stress fluctuations. 

The stress energy tensor of a classical scalar field Q 
with a Lagrangian C and kinetic term 



X 



V„QWQ 



is given by 



It takes the form of a perfect fluid 

T" u = (p + p)U> l U v +p6» v , 
with the association of a 4-velocity pfij 



U„ 



(2JC)V2 



(1) 



(2) 



(3) 



(4) 



which is timelike for the nearly homogeneous cosmolog- 
ical initial conditions. We shall show in WII that it 
remains so throughout structure formation. Here, the 
proper energy density and pressure, or energy density 
and isotropic stress as measured by an observer comov- 
ing with the "fluid" , is 



p = 2XC :X - C, p = C. 



(5) 



Although the mapping from a scalar field to a perfect 
fluid is exact, the fluid equations are not closed without 
a specification of £ to provide the equation of state. 

Scalar fields provide fertile ground for unified dark 
matter models. When combined, the kinetic and poten- 
tial terms possess sufficient degrees of freedom to alter 
the equation of state of the fluctuations separately from 
the background. Both the potential and the kinetic ver- 
sions of UDM posit a separable Lagrangian of the form 



C = L R 



(6) 



where C e and C m behave respectively as dark energy and 
dark matter. Mimicry of a cosmological constant is pos- 
sible if C e is constant or, more generally, slowly varying 
with the field degrees of freedom in comparison to the 
matter piece C m . In what follows, we shall assume that 
the former case holds and set 



PA 



3g 2 » A 
8ttG 



(7) 



The stress energy tensors C e and C m are then sepa- 
rately conserved. Other splittings are of course possible 
(e.g. p e — —C 0) but all choices yield the same observ- 
able predictions. Our choice simplifies the calculation in 
that the dark matter and dark energy pieces interact only 
through gravity. 

The prototype for this sort of unification is an axion- 
likc matter component 



C m = X - m 2 (Q - 



o 2 



(8) 



in a potential whose minimum gives the dark energy den- 
sity C e — — U(Q m i n ). We will refer to this option as 
potential UDM (pUDM). 

Scherrer 19] introduced a similar but alternate ansatz 
based on a modified kinetic term 



Cjri Fi^X A m j n ) , 



(9) 



we call this kinetic UDM (kUDM). The constant F can 
be eliminated in favor of the average matter density and 
kinetic term 



= FX ■ e(4 - 

nun \ 



3e). 



pmO 



3H Q, C 
8irG 



(10) 



where e = (X — A m ; n )/A m j n and "0" denotes evaluation 
at the present epoch. 

The fine tuning between the mass and offset scale in 
pUDM discussed in ^ still exists in kUDM as the coin- 
cidence problem p m o ~ p\. We shall see that kUDM 
also requires eo < 10 -18 (see flVjl . which is similar to the 
pUDM requirement of a purely quadratic potential near 
the minimum or the neglect of self interaction terms in 
the potential. 

Nonetheless, both the pUDM and kUDM model satisfy 
the basic requirement for a successful UDM model in 
that their Lagrangian reaches a non-vanishing minimum 
value C e that can serve as a dark energy component that 
cannot be spatially perturbed. 
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III. EQUATION OF STATE 

Given the mapping between scalar fields and perfect 
fluids described in the previous section, the dynamics of 
the kUDM model are determined by its equation of state, 
the relationship between its pressure and energy density 
in its rest frame. An important property of this purely 
kinetic model is that there is only one degree of freedom, 
the kinetic term X, and therefore the equation of state 
is barotropic to all orders in the density fluctuation, i.e. 
Pm is a function of p m only. 

We begin with the background equation of state 



/ s Pm e(a) 
w m (a) = — = — — ; r 
p m 4 + 3e(a) 



(11) 



To obtain e(a) note that energy conservation provides an 
explicit relation for 



o(e) 



e§(l + eo) 



-1/6 



(12) 



A convenient fitting formula for the required inverse 
relation 



e(a)«e a- 3 U+ (l + e^ 1 ) 



-l\l/3 



-l/ 7 



(13) 



achieves ~ 1% accuracy with 7 = 5/2. Thus, the "mat- 
ter" component behaves as radiation in the background 

1/3 

for a <C e and non-relativistic dark matter thereafter. 

Since the matter component depends only on X, the 
pressure is adiabatic and the fully non-linear equation of 
state 



Pm(pm) = Pn 



/=^e(4 + 3e), 
is barotropic. The adiabatic sound speed 



1 



dp m _ Pm,X _ 1 _ 1 

dp m p m ,x 3 3 VI + / ' 



(14) 



(15) 



then quantifies the pressure response to a density fluc- 
tuation in both the linear and non-linear regimes. Note 
that the pressure is a monotonically increasing function 
of the energy density. Hence, unlike the Chaplygin gas 
model 0, [0| , the matter cannot be destabilized by the 
appearance of high density fluctuations. We will exploit 
this fact in WIl in calculating the abundance of collapsed 
objects. 

Note that the simplification of a purely kinetic UDM 
field is an important feature in the model that is critical 
in the sections that follow. In the more general case (in- 
cluding pUDM), where £ is a function of both Q and X, 
the pressure fluctuations p{X, Q) cannot be expressed as 
a function of p, except in the background. In particular, 




k (Mpc- 



FIG. 1: Modification to the matter transfer function intro- 
duced by kUDM as a function of eo, which quantifies the 
deviation of X from X m i n today. The fluctuations oscillate 
acoustically and decay at small scales in a manner that can 
be conveniently fit by Eqn. 1181 . 



they may be highly time variable. Regardless, the scalar 
field convention is to define the effective sound speed as 
that of the kinetic degree of freedom 0,|2^|, i.e. 



5p 



_ P^X_ 

Ui=0 p t x 



(16) 



since it quantifies the stresses in the frame comoving with 
the field. However, this mapping of field variables onto 
fluid variables is merely formal in the general case. For 
example, the matter component of pUDM acts as colli- 
sionless non-relativistic dark matter even though c~ = 1 
(see Wllfl . In the kUDM model, no such subtlety exists 
and the matter component behaves purely hydrodynam- 
ically. 



IV. LINEAR REGIME 

For the linear fluctuations it is sufficient to note that 
though there is a single kUDM Lagrangian, the matter 
field component is decoupled from the dark energy com- 
ponent and behaves as a fluid. It possesses an adiabatic 
sound speed given by Eqn. (|T5|l with / = 3e(l + 3e/4) and 
negligible anisotropic stress, as the latter scales quadrat- 
ically with spatial field derivatives [see Eqn. (jjJJ]. 

Given that the only gravitational degrees of freedom 
involve the parametrization of the stresses, mimicry of 
ACDM for all linear theory observables above the sound 
horizon, including the cosmic microwave background 
(CMB) and large-scale structure, is guaranteed by con- 
struction p^. Perfect mimicry follows from C e being 
strictly constant, since the spatial fluctuations Sp^ and 
<5pA vanish independently of the scalar field sound speed. 
More generally, allowing C e to depend weakly on Q or X 
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can modify the predictions on the largest scales by intro- 
ducing a relativistic dark energy sound speed [la 1261 127| . 

Below the sound horizon, perturbations in the matter 
component undergo Jeans or acoustic oscillations instead 
of growth, effectively cutting off the linear power spec- 
trum. This behavior places an upper limit on eo such 
that the matter component acts as non-relativistic mat- 
ter during the usual matter dominated epoch a > a oq = 
3 x 10 _4 (f2 m /i 2 /0.14) _1 . Scherrer placed a conservative 
upper limit by demanding that £(0=^ < f or, equiva- 
lent^, e < 4xl0" n (fi m /i 2 /0-14)- 3 0. This constraint 
would ensure that the kUDM behaves as ACDM above 
the horizon at matter radiation equality and thus only 
deep in the linear regime today. To ensure mimicry for 
non-linear structure requires a much more stringent con- 
straint on eo- 

These considerations can be usefully quantified by not- 
ing that the equation of state w m and sound speed drop 
rapidly and monotonically from their relativistic values 
w m = dp m /dp m = 1/3 once a > a(e = 1). After a 
scale exits the sound horizon, matter fluctuations grow as 
CDM independently of scale. Thus, on all larger scales, 
linear fluctuations can be described by a modification to 
the usual transfer function (e.g. pjjj ) 



T(k)=T C vM(k)T Q (k) 



(17) 



To characterize Tq, we solve the coupled Einstein- 
Boltzmann equations that include photo ns, baryons and 
neutrinos in the usual way, following [15l|29|. The result- 
ing transfer function, shown in Fig. ^ is well fit by 



T Q (k) 



with 



V7.74 



3ji(a0 



13 = 0.21 



(-)'' 

\2>AJ 



1/03+1) 



10 



-is 



* L m ,L 

0.14 



0.12 



(18) 



(19) 



Here. 



is the conformal time evaluated at a* = 14e 



1/3 



Note that the power drops by at least an order of mag- 
nitude at x = 5.9, the first extremum of the oscillation. 
For e ^ 0, this drop reflects the lack of the usual loga- 
rithmic growth of the matter fluctuations in the radiation 
dominated era. For larger eo, the fluctuations are further 
suppressed as a function of k due to the change in the 
growth rate in the matter dominated epoch. 

In Fig. |21 we show that the CMB power spectrum pre- 
dictions remain unchanged from ACDM to the cosmic 
variance limit for t < 10 3 and eo < 10~ 16 . The main 
effect of higher e is an enhancement of the radiation 
driving and a reduction of the baryonic modulation of 
the peaks. 



V. HALO ABUNDANCE 

Given that the background pressure of the matter com- 
ponent of kUDM drops as the universe expands, the 




FIG. 2: CMB power spectrum in the fiducial cosmological 
model as a function of eo. For the small eo < 10~ 16 required to 
match early structure, CMB predictions are indistinguishable 
from ACDM. 



largest qualitative difference between this model and 
ACDM in the non-linear regime comes simply from the 
suppression of small scale linear power that is frozen in 
early on. 

This suppression in the initial fluctuations reduces the 
number of low mass dark matter halos in a manner that 
has been well quantified by simulations. For Gaussian 
intial conditions, the number density of halos of mass M 
is controlled by the variance of linear density field pjof 



tT 2 (M;z,e ) 



d 3 k 

^-^P C DM(fc; z)T 2 (k; e )\W(k; M)| 2 . 



(20) 

where W(k; M) is the Fourier transform of the top hat 
window that encloses the mass M at the mean density. 
Specifically, their comoving number density is described 
by the mass function j^J 



dri 



d\nM M 
where v = S c /cr(M) and 



Pmo s dv 



d\nM 



vf(v) = A\j -av 2 [l + [av 2 )~ p ] exp[-az/ 2 /2] 



(21) 



(22) 



A choice of 6 C = 1.69, a = 0.75, p = 0.3, and A such that 
j dvj(y) = 1 fits the results of simulations well. 

At high redshift, where a is small and v large, collapsed 
objects of a given mass are rare. In this regime, the 
number density is exponentially sensitive to c(M; z, eo) 
and hence the amplitude of initial fluctuations. Thus, the 
requirement that there are sufficient low mass halos at 
high redshift to reionize the universe places a constraint 
on kUDM models. 
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FIG. 3: Ratio of mass functions i? mf in kUDM and ACDM. 
The comoving number density of halos of mass M or virial 
temperatures of T v i r is suppressed due to the cut off in the 
power spectrum from linear theory. The suppression increases 
with decreasing mass and increasing redshift. 



Let us define the suppression in kUDM for a given eo 
relative to ACDM in the mass function as 



-Rmf (eo) 



_ dn/dhxM (e ) 
= dn/d\n A/(0) 



(23) 



In Fig. |21 we show this suppression as a function of mass 
for two redshifts representative of reionization and sev- 
eral choices of eo- Note that the suppression increases 
with decreasing mass due to the small-scale cut off in the 
transfer function Tq and with increasing redshift due to 
the exponential sensitivity. 

For the specific question of reionization, the relevant 
quantity for comparison is the fraction of the universe 
that has collapsed into objects of sufficient mass that 
they can cool and form stars. For atomic line cooling, 




FIG. 4: Ratio of collapse fractions i? co i between kUDM 
and ACDM in halos of mass relevant for reionization. Top 
panel: halos that can cool by atomic line transitions T v i r > 
10 4 K. Bottom panel: by molecular hydrogen cooling M > 
10 6 /i _1 Mq. Also shown for reference is the collapse fraction 
itself for ACDM / co i(0). 



the virial temperature of the halo 



T vir = 1.3 x 10 4 K 



(127 



1/3 



M 



l + z 



io / vio 8 ^- 1 ^© 



2/3 



(24) 



must be greater than <~ 10 4 K; for molecular hydrogen 
cooling, M > 10 6 /i _1 Af© but the reionization efficiency 
is lower [3j, |33, HH . The collapse fraction is given by the 
integral over the mass function 



col 



dlnAf- 



M dn 
p m0 dlnM 



(25) 



We will again define the collapse fraction relative to 
ACDM as 



-Rcoi(eo) 



/coi(eo) 



/col(0) ' 

We show this suppression factor in Fig. 



(26) 
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Since even ACDM has difficulty reionizing the uni- 
verse sufficiently early, this ratio cannot be smaller than 
^coi ~ 10 -2 around z ~ 10. Comparing requirements in 
the very similar context of broken scale invariance and 
warm dark matter |32l l33l l35| we place a conservative 
bound of £o 10~ 18 for atomic cooling. If reioniza- 
tion occurs through molecular hydrogen cooling, then the 
limit is even more stringent. Furthermore if indications 
from the CMB polarization measured by WMAP of a 
high reionization redshift z > 10 are confirmed then 
even this regime will be ruled out, given reasonable as- 
sumptions on star formation and the escape fraction of 
ionizing photons |33| . The only loophole is if scale invari- 
ance in the initial power spectrum is strongly broken at 
coincident ally similar scales to compensate. 



VI. HALO SUBSTRUCTURE 

The kUDM model differs from ACDM in the structure 
of dark matter halos both because of the reduction in 
initial small scale power and the fact that it behaves as 
a nearly pressureless fluid instead of a set of collision- 
less particles. These issues will require full cosmological 
simulations to resolve quantitatively, which is beyond the 
scope of this work. However, we shall see that both as- 
pects of the modifications can be reduced to cases that 
have been previously solved in the literature. The latter 
is in qualitative disagreement with some recent observa- 
tions. 

A suppression of small scale power may in fact be desir- 
able to reduce the claimed excess present-day substruc- 
ture in galactic sized halos of the ACDM model. Of 
course, there is tension between these beneficial effects 
and the inability to form enough high redshift halos to 
rcionize the universe. Saturating the conservative bound 
by setting eo ~ 10~ 18 one recovers a suppression of small 
scale power at a similar scale and form to the well-studied 
broken scale invariance model [371 138( . This model was 
explicitly constructed to alleviate the substructure prob- 
lem and so there may remain a small window in kUDM 
parameter space that satisfies these conflicting require- 
ments. Cosmological simulations will be required to ad- 
dress this question quantitatively. Note, however, that 
the substructure problem may have an astrophysical ori- 
gin. 

For the evolution of small scale structure, it is use- 
ful to bear in mind that Eqn. (f 1 4|) implies that the 
pressure remains small compared with the energy den- 
sity until / = 1. Given the low upper bound for eo, 
this corresponds to extremely large overdensities today 
Pmo/Pmo ~ V e o ~ 10 18 . For the relevant / «! case, 
the non-linear equation of state simplifies to 



(27) 



and the sound speed to 
dp r , 



dp r , 



_ 1 Pm 
2 Pm 



(28) 



Note that in this limit the fluctuation to the kinetic term 
is still small, since 



SX 



5p r , 



< 1. 



(29) 



Hence, the fluid 4-velocity [see Eqn. remains time- 
like (see [3!| for a discussion of spacelike kinetic scalar 
fields) and the 3-velocity non-relativistic. Thus, the non- 
linear evolution for kUDM may be calculated with non- 
relativistic hydrodynamic techniques, such as smoothed 
particle hydrodynamics (e.g. [4(1 Elj). 

Qualitatively, density fluctuations will stabilize below 
the Jeans scale, which is the distance sound can travel 
in a dynamical time. The comoving Jeans wavenumber 
becomes 



Y dpm/dpm 



1/2 fVL m h 2 



2\ 1/2 



0.14 



pc- 1 , (30) 



or a physical Jeans scale that is fixed in the matter dom- 
inated epoch for both linear and non-linear density fluc- 
tuations. For a model that has sufficient linear theory 
power to explain high redshift structure and reionization, 
this scale is well below cosmologically interesting scales. 

Given these constraints, the matter component of 
kUDM acts as a pressureless fluid on cosmologically rel- 
evant scales. This type of dark matter has been studied 
extensively in the literature in the context of collisional 
or self-interactin g d ark matter with a vanishingly small 
mean free path [23L |24| . Qualitatively, cold dark mat- 
ter behaves differently from fluid dark matter as soon 
as there exist multiple streams in the dark matter. Since 
cold dark matter is collisionless, it occupies a phase space 
distribution /(x, q, rj), defined in a locally orthonormal 
basis. Its gravitational influence comes through the stress 
energy tensor 



T% = 9 



d 3 q q^q v 
(27T) 3 E(q) 



f: 



(31) 



where q^ is the 4-momentum of the CDM particles 
and g is the degeneracy factor. For CDM in the non- 
linear regime, there can be multiple streams or momen- 
tum states occupied for a given position state. On the 
other hand, kUDM trajectories cannot cross and instead 
kUDM will behave as a fluid and shock. 

Specifically, kUDM or fluid dark matter is challenged 
by observations of a galaxy cluster system observed to be 
undergoing a high velocity merger [42|. Here, the colli- 
sional baryonic gas lags behind the collisionless stars and 
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lensing-observed dark matter due to ram pressure strip- 
ping. Likewise, kUDM is in potential conflict with the 
shape of gravitational wakes of galactic halos moving in a 
cluster halo 0] and the fundamental plane of ellipticals 
in clusters versus the field [4J|. Finally, kUDM would 
not predict as l arg e a triaxiality of halos as measured by 
strong lensing [23, 0> EM Ell 1 although some ambigui- 
ties remain due to projection effects |47J. Note that the 
kUDM is fluid-like even down to scales where the pres- 
sure becomes relevant so that there is no equivalent to 
an adjustable mean free path or cross section to alleviate 
these difficulties |48|. 

VII. DISCUSSION 

Building a successful unified dark matter model is a 
useful exercise, even if its end role is only as a foil 
to ACDM that highlights the successes of the standard 
model. Observations already sharply delineate the basic 
properties of a successful unified model. 

Compared with Chaplygin gas models of unified dark 
matter, kinetic unified dark (kUDM) is more suc- 
cessful as it can explain the well-established observations 
both of the background expansion rate and phenomena 
in the linear regime, such as the cosmic microwave back- 
ground and large-scale structure of the universe. It is not 
without its problems however in the non-linear regime. If 
recent observations in clusters of galaxies are confirmed 
then even this type of model will be ruled out. 

We have shown that the kUDM behaves as a fluid with 
a closed-form barotropic equation of state plus a con- 
stant term that acts as a cosmological constant to all 
orders in structure formation. In linear theory, the high 
density of the matter component in the early universe in- 
troduces a relativistic pressure and hence a cut off to the 
linear power spectrum from Jeans oscillations. We have 
provided a convenient fitting function to the numerical 
results that should assist in future investigations of the 
model. 

In order to produce halos that reionize the universe, 
kUDM must behave like CDM in linear theory on scales 
above k ~ few Mpc . This requirement reduces the vi- 
able kUDM parameter space to the limiting case of highly 



collisional, or self-interacting, non-relativistic dark mat- 
ter with a negligible mean free path [2^, 0] • kUDM is 
then also phenomenologically equivalent to a Chaplygin 
gas UDM, where the pressure is independent of the den- 
sity 01 ( Q ~* 0> see JJ. Such models arc in moderate 
observational conflict with the observed shape of clus- 
ter halos and effects like ram pressure stripping in clus- 
ter substructure. Hydrodynamic simulations of the dark 
matter with the initial conditions provided by the linear 
theory given here may be used to quantitatively address 
these issues in the specific kUDM context. Note that 
these qualitative features are likely to remain valid even 
if the quadratic assumption for the kinetic Lagrangian 
fails far away from the minimum, so long as the implied 
pressure is still of order the energy density once e > 1. 

If these recent observations of clusters of galaxies are 
confirmed only the potential variant of scalar field UDM 
remains viable. Here although the field can still be 
mapped onto a perfect fluid, it possesses a phase space 
structure and mimics CDM down to the de-Broglie wave- 
length of the field, i.e. the kiloparsec regime for m > 
10~ 22 eV ^8|. In this case, the rapidly oscillating equa- 
tion of state makes a direct solution of the fluid equations 
intractable 0, l50| . It can be shown through a WKB 
approximation that the momentum information is car- 
ried by a spatially varying phase to the field, similar to 
a quantum phase space for the wavefunction [5l|. The 
collisionless nature of CDM is then reflected by the lin- 
earity of the wave equation or the quadratic nature of 
the potential. Self-interactions for the pUDM model are 
included by adding anharmonic terms to the potential 

From these examples, it is clear that a successful UDM 
model will likely be as complex phenomenologically as the 
well-tested ACDM model. Complexity reduces the ap- 
peal of the UDM hypothesis, unless a fundamental theory 
can explain the disparate scales introduced by the com- 
ponents that masquerade as the dark matter and dark 
energy. Neither the kUDM or pUDM models in the form 
presented here do so. 
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